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O ■ ABSTRACT 

\m) \ Galaxy-scale gas outflows triggered by active galactic nuclei have been proposed as a key 

fH . physical process to regulate the co-evolution of nuclear black holes and their host galax- 

' ies. The recent detection of a massive gas outflow in one of the most distant quasar, SDSS 

Jl 148+5251 at z = 6.4, presented by Maiolino et al. (2012) strongly supports this idea and 
suggests that strong quasar feedback is already at work at very early times. In a previous 
work, Valiante et al. (2011), we have presented a hierarchical semi-analytical model, GA- 
^ , METE/QSOdust, for the formation and evolution of high-redshift quasars, and we have ap- 

plied it to the quasar SDSS Jl 148+5251, with the aim of investigating the star formation 
history, the nature of the dominant stellar populations and the origin and properties of the 
large dust mass observed in the host galaxy. A robust prediction of the model is that the evo- 
lution of the nuclear black hole and of the host galaxy are tightly coupled by quasar feedback 
in the form of strong galaxy-scale winds. In the present letter, we show that the gas out- 
flow rate predicted by GAMETE/QSOdust is in good agreement with the lower limit of 3500 
M /yr inferred by the observations. According to the model, the observed outflow at z = 6.4 
is dominated by quasar feedback, as the outflow rate has already considerably depleted the 
*^ ■ gas content of the host galaxy, leading to a down-turn in the star formation rate at z < 7 - 8. 

Hence, we predict that supernova explosions give a negligible contribution to the observed 
winds at z — 6.4. 
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1 INTRODUCTION Alatalo et al. 201 1), P-Cygni profiles of FIR OH transitions (e.g. 

Fischer et al. 2010; Sturm et al. 2011), high velocity neutral gas 
The energy released by the black hole (BH) accretion process can . , .. , , . , , . , 

/ , . „ in absorption and high velocity ionized gas (e.g. Nesvadba et al. 

be powerful enough to drive gaseous galactic-scale outflows. Such ,„.„ „ , „ .„ , r c , 

r . . 2010; Rupke & Veilleux 2011) and irom the measurement or ab- 

a negative feedback mechanism is expected to play an important .. . , ... , . . , ~ nnc ,^ ■ , , 

. r r sorption column densities (Arav et al. 2002, 2008) in local quasars 

role in the formation and evolution of quasars (QSOs), by self- .. t im.ii ™no. „■ . i j ■_ j „ T1 

n " J (but see also Chelouche 2008). High velocity and broad [O III] 

regulating the BH growth and eventually quenching star formation, .. , , , ,. ., , . , TT 

° ° c emission observed at z — 2 (e.g. Alexander et al. 2010; Harrison 

therefore affecting the evolution of the physical properties of the . . ~mr,\ i i i ^ a j- 

f „ „ et al. 2012) also trace large-scale gas outflows. Moreover, direct 

host galaxy (e.g. Silk & Rees 1998; Granato et al. 2004; Di Matteo , . . * / „ , ,. ' 

observational evidences ot quasar teedback quenching star iorma- 
etal. 2005; Sprmgeletal. 2005; Ciottietal. 2009; Hopkins & Elvis . . . a , n 

r ° r tion at high redshitt have recently been presented (e.g. Cano-Diaz 

2010; Zubovas & King 2012). Moreover, negative quasar feedback . omo T - . , lnn > 

et al. 2\jl2\ incnas et al. ZU1Z). 

is also required in theoretical models of galaxy evolution to repro- 
duce the space density and properties of old passive and massive Quasar driven outflows have also been observed up to z ~ 6, 
galaxies, observed in the local Universe up to redshift z ~ 2 (e.g. however until recently these were identified only through Broad 
Cimatti et al. 2004; Saracco et al. 2005). Absoiption Lines in their rest-frame UV spectrum (Maiolino et al. 

Observational indications of feedback associated with quasar- 2001, 2004), tracing winds in the vicinity of the accreting black 
driven massive outflows come from the detection of broad wings of hol e and accounting only for a tiny fraction of the amount of gas in 
molecular (CO and HCN) emission lines (e.g. Feruglio et al. 2010; the whole host galaxy. 

In a previous work (Valiante et al. 2011, hereafter VI 1) we 
have presented a semi-analytical model, GAMETE/QSOdust, for 
* E-mail: valiante@arcetri.astro.it the formation and evolution of high redshift quasars. This model 
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Table 1. Observed and inferred physical properties of the z = 6.4 QSO Jl 148 and its host galaxy. See Valiante et al. 201 1 for details. 



Quantity 


Description 


Value 


Reference 


M 1450A 


Continuum rest frame AB magnitude 


-27.82 


Fan et al. 2003 


M h 


Mass of the host dark matter halo 


10 13 M© 


Fanet al. (2004); VI 1 


Mbh 


Mass of the super massive black hole 


3+ 3 jj X 10 9 M Q 


Barth et al. 2003; Willott et al. 2003 


Lfir 


Far Infrared Luminosity 


(2.2±0.3)L G 


VI 1 and references therein 


M dyn 


Mass of the dynamical gas for a inclination angle i = 65° 


(5.5 ±2.75) X 10 10 M 


Walter et al. 2003 


M H2 


Mass of the molecular gas from the CO observations 


1.6 x 10 10 M Q 


Walter et al. 2003 


M star 


Stellar mass computed as Mj v „ - Mjj 2 


(3.9 ±2.75) x 10 10 M Q 


Walter et al.2003 


M star 


Stellar mass from M BH /M 5lar ~ 0.002 


~ 2.14 X 10 12 M 


Marconi & Hunt 2003 


SFR 


Star formation rate 


(3.8 ± 0.57) X 10 3 M Q /yr 


Bertoldi et al. 2003 


Z 


metallicity from the observations of the narrow line regions in QSOs 


1.32 , 1Q Ad 


Matsuoka et al. 2009 


M lius , 


Mass of dust 


3.4+[- 3 * X 10 s M 


VI 1 and references therein 


M out fi 


quasar-driven gas outflow rate 


> 3500 M /yr 


Maiolino et al. 2012° 



" The gas outflow rate has been inferred from the observations of Jl 148 only after the model presented in VI 1 was published (see text) 



has enabled us to constrain the star formation history (SFH) and 
the nature of the stellar populations (e.g. the stellar initial mass 
function, IMF) of the well studied high redshift quasar SDSS 
Jl 148+5251 (hereafter Jl 148) at z=6.4, by linking the evolution 
of the nuclear black hole (BH) with the chemical properties of the 
host galaxy. The main observed and inferred physical properties of 
J 1148 are summarized in Table [T] In VI 1 we have found that the 
evolutionary scenarios that allow us to reproduce the physical prop- 
erties of the quasar J 1148 and of its host galaxy, such as the mass 
of the nuclear BH, the dynamical, gas and dust masses, are tightly 
constrained. More specifically, if stars are assumed to form accord- 
ing to a standard IMF, such as a Larson IMF with a characteristic 
mass of m c h = 0.35 M Q , the total stellar mass exceeds the upper 
limit set by the observed dynamical mass by a factor of 3-10, de- 
pending on the adopted star formation model (see VI 1 for further 
details). Alternatively, the total stellar mass can be accommodated 
within the dynamical mass limit if a top-heavy IMF is assumed for 
all stellar populations formed in the host galaxy, with a character- 
istic stellar mass of m c h = 5 M Q . 

More important, for the purpose of the present investigation, a 
fundamental prediction of the model was that "a powerful outflow 
is launched by the quasar during the latest ~ (100 - 200) Myr of 
the evolution". This is a natural consequence of continuous neg- 
ative feedback that operates along the merger history of the host 
galaxy and it is required to regulate BH growth, star formation and 
to reproduce the observed dynamical and chemical properties in all 
of the evolutionary scenarios presented in VI 1. 

Recently, Maiolino et al. (2012) have detected, for the first 
time, a massive gas outflow, ascribed to strong quasar feedback, 
through IRAM PdBI observations of the [CII] 158yum transition 
in the host galaxy of J 1148. They have revealed broad extended 
[CII] wings which indicate a large mass (7 x 10 9 M Q ) of outflowing 
atomic gasAssuming a spherical outflow and that maximum veloc- 
ity observed in the wings, v = 1300km/s is the de-projected veloc- 
ity of the outflow, they obtain a conservative lower limit on the out- 
flow rate of 3500 M G /yr, which represents the highest outflow rate 
ever detected (see Maiolino et al. 2012 for further details on the ob- 
servations). The associated kinetic power is barely consistent with 
the kinetic power that can be produced by a starburst driven wind, 
even assuming a star formation rate of ~ 3000M o /yr as inferred 
from the strong far-IR thermal emission detected by (sub-)mm ob- 
servations of J1148 (Bertoldi et al. 2003; Beelen et al. 2006). On 
the other hand, the kinetic energy is about 0.6% of the bolomet- 
ric luminosity of the quasar, meaning that the observed outflow is 
more likely to be powered by the active galactic nucleus (AGN). 

Here we show that this observational result is a strong confir- 



mation of the theoretical models presented in VI 1. In fact, in all of 
the evolutionary scenarios that successfully reproduce the observed 
properties of Jl 148, the outflow rate at z = 6.4 is largely dominated 
by the AGN and it is consistent with the observed lower limit. 

In what follows, we adopt a Lambda Cold Dark Matter 
(ACDM) cosmology with £l m = 0.24, Cl A = 0.76, Cl b = 0.04, and 
Ho = 73 km/s/Mpc. The age of the Universe at a redshift z = 6.4 is 
900 Myr. 



2 MODEL DESCRIPTION 

Through a binary Monte Carlo algorithm with mass accretion based 
on the EPS theory, we have simulated 50 hierarchical merger histo- 
ries of a 10 13 M o dark matter (DM) halo at z = 6.4. As discussed in 
VI 1, a [10 12 - 10 13 ] M G host DM halo is required to match the ob- 
served space density of z ~ 6 QSOs. The dependence of the results 
on the DM halo mass will be discussed in section[5] Using the code 
GAMETE/QSOdust, we have followed in a self-consistent way the 
build-up of the Jl 148 central SMBH and of its host galaxy along 
these formation paths (VI 1). 

The nuclear SMBH forms starting from 10 4 /r' M Q seed BHs 
hosted in very high redshifts halos. Along each merger tree, the 
assembly of the host galaxy and its SFH are driven by binary merg- 
ers and mass accretion. Similarly, seed BHs grow by coalescences 
with other BHs and by gas accretion. The co-evolution of these two 
components, the SMBH and its host galaxy, is mostly controlled by 
quasar feedback in the form of a galactic-scale wind, which self- 
regulate the BH growth and eventually halts star formation. 

The evolution of the mass of gas, stars and metals in the inter- 
stellar medium (ISM) is followed by taking into account infall and 
outflow processes from/to the external medium. Metal and dust en- 
richment is implemented in all progenitor galaxies by considering 
the contribution from stars of different masses, such as Asymptotic 
Giant Branch (AGB) stars and Supernovae (SNe), taking into ac- 
count the proper stellar lifetimes. Moreover, subsequent dust grain 
re-processing in the ISM, namely destruction by interstellar shocks 
and grain growth in molecular clouds (MCs), is also taken into ac- 
count. 

Finally, different SFHs for the QSO host galaxy are explored: 
quiescent models, where the efficiency of star formation is indepen- 
dent of galaxy mergers, and bursted models in which the star for- 
mation efficiency is enhanced during major mergers between pro- 
genitor galaxies. 
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3 MODEL RESULTS 

In this section we present the models that in VI 1 have been identi- 
fied to reproduce the available observational constraints on J 1 148. 
In all models, the free parameters, namely the efficiency of star for- 
mation, BH accretion and quasar-driven wind have been chosen to 
reproduce the BH and gas mass of J 1 148. 

As discussed in VI 1, the inferred stellar mass, a crucial in- 
formation for the evolutionary models, is affected by large uncer- 
tainties. This is the reason why two different values are given in 
Table Q] and the models that we have explored adopt progressively 
higher star formation (SF) efficiencies, and therefore larger final 
stellar masses. We thus, investigate the corresponding effects on 
the evolution of the chemical properties of host galaxy (metal and 
dust content). 

Four out of the several models presented in VI 1 have been 
found to successfully reproduce the properties of the QSO Jl 148: 
Ql t.h., Bl t.h., Q2 and B3. The first two are quiescent (Ql t.h.) 
and bursted (Bl t.h.) star formation models where stars form ac- 
cording to a top-heavy IMF, while the quiescent (Q2) and bursted 
(B3) models assume a standard IMF 




10 7 10 8 10 9 10 10 10 11 
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10 12 1Q 13 



3.1 BH-M stor relation 

For convenience, we report in figure [T] the evolution of the 
M BH -M s , ar relation predicted by the four models introduced above. 
The location of II 148 on the M BH -M s , ar plane is indicated by the 
solid circle (with M !fo , given by the observations, see Table[T) and 
is compared with the M BH -M slclr relation observed in local quasars 
and galaxies (open squares) and the empirical fit taken from Mar- 
coni & Hunt (2003). A detailed discussion can be found in VI 1 . By 
construction, all models reproduce the mass of the SMBH in II 148 
at z=6.4, but predict different final stellar masses going from the 
value inferred from the observations, ~ 3.9 x 10 10 M Q (models Ql 
t.h. and B 1 t.h.) up to values progressively closer to the observed 
local relation, ~ (1 - 4) x 10" M Q (Q2 and B3, respectively). 



3.2 Chemical properties and SFH 

By using the chemical evolutionary model with dust described in 
VI 1, we have analyzed the dependence of the predicted Jl 148 
chemical properties (in particular of metals and dust) on both the 
SFH (quiescent vs bursted star formation rates, SFRs) and IMF 
(standard vs top-heavy). The results of the four selected models 
are shown in figure [2] In this plot, we present the predicted evolu- 
tion of the mass of gas, stars, metals and dust compared with the 
available observations for Jl 148. As it can be seen from the figure, 
independently of the SFH and IMF, the mass of gas and metals both 
grow up to a maximum value, reached at redshift z ~ 8, and than 
decrease due to the strong effect of the quasar feedback. 

The SFR of high redshift quasars is usually estimated from the 
FIR luminosity, adopting the L Fm - SFR relation derived by Ken- 
nicutt (1998) under the assumption that the dominant dust-heating 
mechanism is radiation from young stars. The conversion factor 
usually adopted to convert L FjR to SFR, ~ 5.8 x 10 9 L Q , has been 
derived assuming a 10 - 100 Myr old burst of star formation and 
a Salpeter IMF (see Kennicutt 1998). Following this prescription, 
the FIR luminosity estimated for 11148 (see Table Q} gives a SFR 
~ (3.8 + 0.57) x 10 3 M /yr (Bertoldi et al. 2003). However, it has 
been pointed out that the inferred SFR would be about an order of 
magnitude lower (~ 180 M Q /yr) if the Schmidt-Kennicutt law is 
instead adopted (Dwek et al 2007; Li et al. 2008; VI 1). 



Figure 1. The evolution of the Jl 148 BH mass as a function of the stel- 
lar mass for quiescent (upper panels, Ql t.h. with black line and Q2 with 
azure line) and bursted models (lower panel, B l t.h. gold line and B3 blue 
line). The filled circle and open squares are the Jl 148 and local QSOs and 
galaxies observed Mbh-M jmi - relations, respectively. The dashed line is the 
Marconi & Hunt empirical fit to the local relation. Solid lines are the four 
different models that in V 1 1 have been selected to reproduce the available 
observations of Jl 148 (see text), averaged over 50 different merger tree re- 
alizations. The shaded areas represent the 1 cr dispersion. 
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Figure 2. The ISM chemical evolution of the Jl 148 host galaxy in the four 
different models which best reproduce the observed properties. Solid lines 
are the evolution of the mass of gas (blue), stars (yellow), metals (green) 
and dust (red) averaged over 50 hierarchical merger histories with the lcr 
dispersion given by the corresponding shaded areas. The blue arrow indi- 
cates the molecular gas mass (M# 2 ), considered as a lower limit for the total 
gas mass, the yellow solid square is the stellar mass computed as the differ- 
ence between the dynamical and molecular gas mass, the red triangle is the 
mass of dust and the green region represents the mass of metals (Z X Mh 2 )■ 
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In this work we have further revised the SFR inferred from the 
observations, by taking into account the different IMFs adopted in 
our models and we rescale the L F1R - SFR conversion factor for 
a Larson IMF with m ch = 0.35 M Q (standard IMF models) and 5 
M (top-heavy IMF). With this correction we obtain star formation 
rates ~ 3 to 30 times lower, more specifically: SFR ~ (1087 ± 
163)M G /yr, for a standard IMF, and SFR ~ (1 13 ± 17) M Q /yr, for a 
top-heavy IMF. These values are still consistent with results quoted 
in previous works (Maiolino et al. 2005; Dwek et al. 2007; Li et 
al. 2008). Yet, as discussed in VI 1, we use these values as an upper 
limit to the effective rate of star formation in Jl 148, since the active 
quasar itself may give a non negligible contribution to dust heating 
(Bianchi, Valiante & Schneider in preparation). The quiescent and 
bursted SFHs for the four selected models are shown in the left 
panels of figures [3]and[4] 



4 QUASAR AND SN-DRIVEN GAS OUTFLOW 

In GAMETE/QSOdust, galactic outflows can be driven by SN ex- 
plosions and by the AGN. Both stellar and quasar feedback are con- 
sistently modelled in the form of energy-driven winds sweeping the 
surrounding material away from the galaxy. The total gas outflow 
rate is given by dM eJ (t)/dt = dM e;>SA i(r)/df + dM e ^ aN {t)ldt, with: 



Models with a standard IMF (Larson, m ch -0.35 



dt 



= 2e w ,AGN£r( — ) Mic 

Ve 



and 
dM, 



ej.SN 



2e w E s 



dt 



-*sn(0, 



(1) 



(2) 



where is the escape velocity of the galaxy, Rsn(1) is the SN ex- 
plosion rate, M accr is the BH gas accretion rate, e r is the radiative ef- 
ficiency, fixed to be 0. 1 and £sn is the average supernova explosion 
energy, assumed to be 2.7 x 10 52 erg for PISNe and 1.2 x 10 5 ' erg 
for core collapse SNe. The efficiency of SN-winds (e w = 2 x 10" 3 ) 
has been calibrated by the chemical evolution model applied to the 
Milky Way and its dwarf satellites (Salvadori et al. 2007, 2008). 
The AGN- wind efficiency {e Vi AGN = 5 x 10~ 3 ) and the BH accre- 
tion efficiency (a = [ 1 80 - 200] , according to the particular model), 
have been instead chosen to reproduce at the same time the gas and 
BH mass of J1148. 

The right panels of figures |3]and|4]show the redshift evolution 
of the gas outflow rate predicted by the selected models, as labelled 
in the figures, by showing the separate contribution of SNe and the 
quasar. As it can be seen from these figures, for z < 12 the QSO 
outflow rate exceeds by more than two orders of magnitude the SN 
outflow rate in all the different models. 

For all models, the predicted outflow rate at z = 6.4 is fully 
consistent with the lower limit to the gas outflow rate inferred 
by Maiolino et al. (2012). With this large outflow rate, the host 
galaxy should be rapidly cleared of its gas content and star forma- 
tion should be quenched. Indeed, the predicted SFRs systematically 
show a down-tum at redshifts z < 7 - 8, when the effect of the 
quasar-driven wind on the evolution of the ISM gas mass becomes 
dramatic (see blue lines in figures |2}. Such a decline is reflected 
also into the SN rate and, therefore, in the SN-driven outflow rate, 
which is orders of magnitude lower than the quasar-driven one. 

The formulation adopted in VI 1 for the BH accretion rate en- 
sures that the quasar feedback is active over almost the entire life- 
time of the host halo. This feedback results into efficient gas ejec- 
tion at z < 8, when the BH mass and accretion rate are higher (see 
VI 1 for a detailed discussion). 
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Figure 3. The star formation (left panels) and the gas outflow (right pan- 
els) rates as a function of redshift, predicted by models B3 and Q2. In both 
models a Larson IMF with a characteristic mass m c i, = 0.35 M Q has been 
adopted. All the curves are the averages over 50 random merger tree real- 
izations of the quasar host galaxy, with shaded areas representing the lcr 
dispersion. Blue dashed and red solid lines in right panels represent the 
mass of gas ejected per unit time by SN and quasar driven winds, respec- 
tively. The arrows indicate the upper limit to the SFR of J 1 148, corrected 
for our adopted IMFs or to the gas outflow rate, obtained by Maiolino et al. 
(2012). 
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Figure 4. The same as in figuref5]but for models Bl and Ql. In these two 
models a Larson IMF with m c i, = 5 M G has been adopted. 
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5 CONCLUSIONS 

In this work, we have shown that the formation history of SMBH 
at the center of high redshift quasars is strongly linked to the prop- 
erties of their host galaxies. 

We use a state-of-the-art model (VI 1) that allows us to simu- 
late a large number of independent hierarchical histories of quasar 
host galaxies, following at the same time the star formation histo- 
ries, chemical evolution, and the gradual build-up of the nuclear 
black hole by mergers and gas accretion. 

We have recently applied this numerical model to one of the 
most distant quasars, SDSS 11148 at z = 6.4 and we have shown 
that the joint constraints provided by the dynamical, gas, and dust 
masses together with the mass of the nuclear SMBH, allow to re- 
construct the star formation histories and the nature of the dominant 
stellar populations in these systems. A robust prediction of these 
models is that the evolution of the nuclear black hole and of the 
host galaxy are tightly coupled by quasar feedback in the form of 
strong galaxy-scale winds. 

The recent detection of broad wings in the [CII]158yum line 
associated to SDSS 11148 by Maiolino et al. (2012), is a strong 
evidence of clear quasar feedback in the early Universe. 

The comparison between model predictions and observational 
data shows that the gas outflow rates predicted by the models are 
fully consistent with the lower limit (3500 M Q /yr) inferred from the 
observations. In addition, it lends strong support to the interpreta- 
tion that the dominant driving mechanism is quasar feedback. Ac- 
cording to the physical prescriptions and parameter values adopted 
in the models presented here, the outflow rate has already signifi- 
cantly affected the gas content of the galaxy, leading to a down-turn 
in the star formation rate at z < 7 - 8. 

It is important to stress that the theoretical findings have been 
derived independently, well before the observations. Therefore, 
strong quasar feedback has been a natural by-product of models 
that have been identified to reproduce the observed dynamical and 
chemical properties of the host galaxy. In all models, the SN rate 
is constrained by the SFR inferred by the observations and by the 
chemical properties of the host galaxy. The associated SN-wind 
efficiency has been calibrated using the same chemical evolution 
model to reproduce the global properties of the Milky Way Galax- 
ies (Salvadori et al. 2008). To be at least comparable with the ob- 
served outflow rate, the SN wind efficiency should be higher by 
~ 3 - 4 orders of magnitude. There are no physical motivations to 
justify such an enormous efficiency. 

All the proposed scenarios, i.e. different combinations of 
SFHs and IMFs, provide consistent results for the quasar-driven gas 
ejection rate. Thus, the recent observation reported by Maiolino et 
al. (2012) does not help to break the degeneracy between the pro- 
posed models: accurate observations of the stellar mass are nec- 
essary to further constrain the properties of the Jl 148 host galaxy 
therefore helping in the identification of the most appropriate IMF 
and SFH. 

Finally, the comparison of our predicted gas outflow rate with 
that observed for II 148 by Maiolino et al. (2012) allow us to put ad- 
ditional constraints on the mass of the DM halo hosting the SMBH. 
With a DM halo mass equal to 10 12 M o , both the BH accretion and 
feedback efficiencies need to be reduced in order to reproduce the 
observed BH and gas masses at z = 6.4. In fact, a smaller DM halo 
mass implies a reduced gas mass available at each redshift along 
the hierarchical history. As a result, the final outflow rate is more 
than one order of magnitude smaller than for the reference model 



with DM halo mass of 10 13 M o and inconsistent with the observed 
lower limit. 
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